Functional renormalization group approach to the Ising-nematic quantum critical 

point of two-dimensional metals 



o 

(N 






I 

-US 

o 
o 



> 

in 

rn 

O 

> 



Casper Drukier, Lorenz Bartosch, Aldo Isidori, and Peter Kopietz 

Institut fur Theoretische Physik, Universitdt Frankfurt, 

Max-von-Laue Strasse 1, 60438 Frankfurt, Germany 

(Dated: March 12, 2012) 

Using functional renormalization group methods, we study an effective low-energy model describ- 
ing the Ising-nematic quantum critical point in two-dimensional metals. We treat both gapless 
fermionic and bosonic degrees of freedom on equal footing and explicitly calculate the momentum 
and frequency dependent effective interaction between the fermions mediated by the bosonic fluctu- 
ations. Following earlier work by S.-S. Lee for a one-patch model, Metlitski and Sachdev [Phys. Rev. 
B 82, 075127] recently found within a field-theoretical approach that certain three-loop diagrams 
strongly modify the one-loop results, and that the conventional 1/A'^ expansion breaks down in this 
problem. We show that the singular three-loop diagrams considered by Metlitski and Sachdev are 
included in a rather simple truncation of the functional renormalization group flow equations for 
this model involving only irreducible vertices with two and three external legs. Our approximate 
solution of these flow equations explicitly yields the vertex corrections of this problem and allows 
us to calculate the anomalous dimension rj^p of the fermion field. 

PACS numbers: 05.30.Rt, 71-lO.Hf, 71.27.-fa 



I. INTRODUCTION 

Inspired by the puzzling normal-state properties of 
the copper-oxide superconductors, the search for possi- 
ble non-Fermi hquid states of metals continues to be a 
central topic in the theory of strongly correlated elec- 
trons. An interesting new clue toward an understanding 
of the cuprates and other materials comes from exper- 
iments on a variety of materials, indicating a nematic 
phase transition.'2'Zl This is a quantum phase transition 
which, while preserving translational symmetry, breaks 
the lattice rotational symmetry from square to rectan- 
gular, i.e., the invariance under rotations of the system 
in the x-y plane by 90° is lost. At the quantum criti- 
cal point, electrons couple strongly to order parameter 
fluctuations, leading to a destruction of the Fermi liquid 
state. The resulting distortion of the Fer mi su rface is 
also referred to as a Pomeranchuk transition.lSEl 

The conventional theoretical approach to quantum 
critical phenomena is the so-called Hertz-Millis approach, 
where the interaction between electrons is decoupled via 
a (bosonic) Hubbard-Stratonovich transformation and 
the fermionic degrees of freedom are integrated out.'iSHini 
However, as in the metallic state the electrons are gapless, 
this approach usually leads to singular vertices, which es- 
pecially in low dimensions need to be treated with care. 
It can therefore be advantageous not to integrate out the 
fermions at all. 

As concerns the nematic phase transition, the most ef- 
fective scattering processes of electrons take place when 
the momentum of the bosons is locally almost tangential 
to the Fermi surface. The phase transition can therefore 
be modeled by coupling electrons in the vicinity of two 
patches of a Fermi surface to a gapless scalar Ising order 
parameter field. Indeed, the coupling of gapless fermions 
to gapless bosonic fluctuations is known to give rise to 



non-Fermi liquid behavior.li^E^l The corresponding field 
theory is very similar to a nonrelativistic gauge theory, 
which has been studied intensively,'!^"!^! starting with the 
important work by Holstein, Norton, and Pincus.^^ Such 
gauge theories have applications in a number of different 
problems, such as the description of the half- filled Landau 
level,^ the description of spin liquids in terms of spinous 
which can form a critical spinon Fermi surface, ^^ '^^ or the 
description of the instability of a ferromagnetic quantum 
critical point.l^il Similar gauge theories have also been 
used to describe fermions on a honeycomb lattice inter- 
acting through an electromagnetic gauge field. ^^ ^^ In all 
cases, the low-energy behavior is expected to be described 
by a scale-invariant scaling theory. The si ngle-particle 
Green function G{uj, k) was calculatecP^I^Sl for a spher- 
ical Fermi surface within the random phase approxima- 
tion (RPA), resulting in 
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Here, uj and k are the frequency and momentum of the 
electron, A^ = A'sgn {cu) + iA" is a complex constant 
depending on the sign of w, and ^k = ^F(|fc| — ^f), where 
vp is the Fermi velocity and kp is the Fermi momentum, 
denotes the single-particle excitation energy. While the 
static part of the self-energy remains unrenormalized, its 
dynamic part implies that both the renormalized energy 
and damping rate of the electron scale in exactly the 
same way. Consequently, there are no well-defined sharp 
quasiparticles and Landau's Fermi liquid theory breaks 
down. 

For a long time, it was thought that the above scenario 
holds true when going beyond the RPA.^^ It was believed 
that this can be justified by considering the limit of large 
N, where N is the number o f ferm ion flavors. However, 
it was recently shown by Lee^^Sll that, even if one con- 
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FIG. 1. These three- loop diagrams have been identified by 
Metlitski and Sachdev^^'^ to give singular corrections to the 
one-loop results for the bosonic and fermionic self-energies. 
Solid arrows denote the fermionic single-particle Green func- 
tions within the one-loop approximation, while wavy lines rep- 
resent the RPA propagator of the bosonic fluctuations. The 
black dot is the bare interaction vertex between one boson 
and two fermion fields. The diagrams in (a) are three-loop 
corrections to the fermionic self-energy, while those in (b) rep- 
resent corrections to the bosonic self-energy of the so-called 
Aslamazov-Larkin type. 



siders only scattering processes in the vicinity of a sin- 
gle patch of the Fermi surface, the coupling to a gapless 
gauge field results in very strong correlations such that, 
even in the large-A'^ limit, the theory remain s stro ngly 
coupled. Subsequently, Metlitski and Sachdev^^ffl con- 
sidered a more realistic two-patch model where the two 
patches of the Fermi surface to which a given bosonic 
momentum is tangent are retained; they derived a scal- 
ing theory and explicitly calculated corrections to the 
bosonic and fermionic self-energies up to three loops, us- 
ing the one-loop propagator in internal loop integrations. 
Metlitski and Sachdev identified certain three-loop con- 
tributions to the bosonic and fermionic self-energies (see 
Fig. II]) , which give rise to logarithmically divergent cor- 
rections to the one-loop RPA results. Exponentiating 
these logarithmic terms, they obtained the following ex- 
pression for the retarded propagator of the fermions. 
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where ry^ is the anomalous dimension of the fermion 
field and z is the fermionic dynamic critical exponent. 
Using scaling relations for the fermionic and bosonic 
Green functions, Metlitski and Sachdev argued that the 
fermionic dynamic exponent is given by z = Zb/2, where 
Zb is the corresponding bosonic dynamic exponent. Their 
explicit calculations^ show that z;, = 3 is not renor- 
malized by fluctuations up to three loops, implying that 
z = 3/2 is correctly given by the one-loop approxima- 
tion. On the other hand, for the fermionic anomalous 
dimension, Metlitski and Sachdev obtain the finite re- 
sult rjjf, = 0.068 at the Ising-nematic transition, whereas 
r]^ = within the one-loop approximation. 

On a technical level, the reason for the breakdown of 
the large-A^ expansion can be traced back to the fact 
that the curvature of the fermion propagator comes with 
a factor of N. This in combination with a cancellation 
of the curvature in a set of planar diagrams eventually 



leads to the breakdown of the large- A^ expansion. In the 
words of Chubukov,!^ there is hidden one- dimensionality 
in two-dimensional systems. 

Albeit 1/iV can not be used as a control parameter, it 
was suggested by Mross et al. ^ following earlier work at 
finite N by Nayak and Wilczek,'321 to use z^ as a tunable 
parameter. In this case it is possible to consider the lim- 
its iV —> oo and Zfe — 2 — > while keeping the product 
N{zb — 2) finite to bring the calculation under control. 
Using a tunable z^ is a sensible strategy because a non- 
local interaction is not expected to be renormalized. Ex- 
trapolating the results obtained by Mross et al^^ to the 
physically relevant case Zf, = 3 and N = 2, one obtains 
for the anomalous dimension of the fermion field rj^, « 0.6 
(using our definition (1.2) of 77^). Obviously, this value 
is much larger than the estimate ry^ ~ 0.068 by Metlitski 
and Sachdev.!^ 

Even though the calculations in Refs. UHl and [32] are 
based on the field-theoretical renormalization group, the 
fact that two independent calculations involving differ- 
ent types of approximations produce different values for 
T]^ shows that on a quantitative level there are still open 
questions. Due to the sign problem in quantum Monte 
Carlo calculations and the fact that dynamical mean field 
theory can essentially only predict mean field exponents, 
the number of alternative methods to verify the correct- 
ness of the anomalous scaling properties of the Ising- 
nematic transition is limited. In this work, we study 
this problem by means of the one-particle irreducible 
implementation of the functional renormalization group 
(FRG) method,'231t23 -which is a modern implementation 
of the Wilsonian renormalization group idea. The fiex- 
ibility of FRG methods to deal with systems involving 
both fermionic and bosonic fields has already been used 
by several authors .'^SHlll In particular, in Refs. Hnms", and 
HH it has been shown that it can be advantageous to intro- 
duce a cutoff parameter A which regularizes the infrared 
divergences only in the momentum carried by the bosonic 
field (the momentum-transfer cutoff scheme). We show 
in this work that this cutoff scheme is also convenient to 
study the nematic quantum critical point. 

The rest of this work is organized as follows. After 
introducing the model system and defining our notation 
in SecUD we give in Sec. |III| the FRG flow equations for 
the self-energies and vertex corrections in general form. 
We also introduce the momentum-transfer cutoff scheme 
and show that the singular three-loop diagrams shown 
in Fig. [1] are contained in a rather simple truncation of 
the hierarchy of FRG flow equations involving only irre- 
ducible two-point and three-point vertices. In Sec. |IV[ 
we show how to recover the known one-loop results for 
the momentum- and frequency-dependent fermionic and 
bosonic self-energies by integrating the FRG flow equa- 
tions ignoring vertex corrections. In the main part of 
this work, given in Sec. Ivj we consider the system of 
FRG flow equations including vertex corrections. We ex- 
plicitly calculate the effect of vertex corrections on the 
value of the fermionic anomalous dimension, rj^, to lead- 



ing order in the small parameter Zb — 2. In the con- 



cluding section VI , we summarize our results and discuss 
some open problems. We have added two appendices 
with more technical details. In Appendix A we derive 
skeleton equations relating the purely bosonic two-point 
and three-point functions to the fermionic propagators 
and irreducible vertices. These skeleton equations are 
useful to close the infinite hierarchy of FRG fiow equa- 
tions. Finally, in Appendix B, we explicitly evaluate the 
symmetrized fermionic loop with three external bosonic 
legs (the symmetrized three-loop vertex) for our model 
system. 



II. DEFINITION OF THE MODEL 

We are interested in a minimal model describing 
the coupling of electrons in the proximity of a two- 
dimensional Fermi surface to a gapless scalar Bose field. 
In particular, the bosonic field can describe the fiuc- 
tuations of a scalar order parameter near the onset of 
a metallic Ising-nematic phase, such as a d-wave ne- 
matic state in a two-dimensional square lattice, where 
the point-group symmetry of the lattice is reduced from 
square to rectangular. However, the gapless scalar Bose 
field can also describe the fluctuations of an emergent 
gauge field minimally coupled to a two-dimensional Fermi 
surface. For example, this can be physically realized in a 
system where a spin-liquid phase is described in terms of 
fermionic degrees of freedom (spinous), thereby causing 
the emergence of a U{1) gauge symmetry in the system: 
The critical fluctuations near the onset of the spin-liquid 
phase are then described by the coupling of the spinon 
Fermi surface to the corresponding f/(l) gauge fleld. The 
general form of the action for our model can be written 



asS = S., 
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where '0 and cf) denote two-dimensional Fermi and Bose 
fields, respectively, and O[0,'0] is a bilinear operator 
in the fermion fields which has the same symmetry as 
the order parameter field (j). In Eqs. (2.1)-(2.3), K — 



(iuj, k) denotes fermionic Matsubara frequency and two- 
dimensional momentum, while K = (zw, k) denotes the 
corresponding bosonic quantities. The integration sym- 
bols are defined by /^ = {PV)~^ J^w Sfcj ^^'^ similarly 
for the bosonic quantities, where /3 is the inverse tem- 
perature and V is the volume. Throughout this work 
it is understood that we eventually take the zero tem- 
perature limit (/3 — >■ oo) and the infinite volume limit 
(V — )■ oo). The index a = 1, . . . , A^ labels N different 
flavors of the fermion field. The fermionic energy dis- 



persion ^fc is defined relative to the Fermi energy ep, i.e. 
^k = ^k—^F, while, in the bosonic dispersion, po plays the 
role of a mass (or gap) term which measures the distance 
to the quantum critical point: At the quantum critical 
point, Po — 0, such that order parameter fluctuations 
become gapless. The absence of higher-order terms in (p 
and gradients of in the action defining our model can 
be justified by a dimensional analysis, which shows that 
such higher-order terms become irrelevant at the critical 
point. 



In the general form given in Eqs. ( 2. 1 )-( 2.3 ) , the action 



of our model is still too complicated to be treated ana- 
lytically with renormalization group or field-theoretical 
methods. However, as pointed out by Metlitski and 
Sachdev22 (gge also Ref. ^35), the relevant critical fluc- 
tuations can be described by a simplifled minimal action 
involving only fermion fields with momenta close to two 
opposite patches on the Fermi surface. The reason is 
that the most singular scattering processes mediated by 
a given bosonic mode with momentum k involve only 
fermions lying on patches of the Fermi surface which are 
almost tangential to the bosonic momentum k. The sit- 
uation is shown graphically in Fig. [2J where the label 
a = ±1 denotes the two patches of the Fermi surface 
which are tangential to a given bosonic mode with mo- 
mentum parallel to fcj^ in the figure. In order to describe 
the singular behavior of the fermionic and bosonic Green 
functions at the critical point we can therefore restrict 
the general model involving fermions on the whole Fermi 
surface to a so-called two-patch model, characterized by 
the following Euclidean action: 

Vtches[Vi,V',0] = sSM + s^ii^ + Sii^AA], (2.4) 



SSM - - / Y.^^^ ~ ^kWK.rKa, (2.5) 

So[4'] = lJ^f^'<l^-K^K. (2.6) 

JKi JK2 JKs „_^ 

xrf'^°nifi;i^2;i^3)v;^,>K..'/'K3-(2.7) 



In the above expressions the fermion fields are now char- 
acterized by an additional upper index a — ±1 labeling 
the two patches on the Fermi surface centered at the two 
opposite momenta kp = akp, as shown in Fig. [2] By 
construction, the momenta of the fields characterizing 
the two-patch model are intended to lie in the vicinity 
of the Fermi momenta kp, so that |fc||| and |fcj^| (the 
momenta relative to the Fermi momenta, as shown in 
Fig. I2| should be much smaller than \kp\. In particular, 
one should impose a cutoff, A_l ~ kp/^9 on the momenta 
perpendicular to the Fermi surface normal, where A0 is 
the angular extension of the patch. However, as long as 
integrals over such momenta turn out to be ultraviolet 
convergent, we can effectively send this cutoff to infinity 




FIG. 2. The two-patch model considered in this work involves 
only two types of fermion fields with momenta close to two 
opposite patches on the Fermi surface centered at itkr- The 
fermionic momenta are measured locally with respect to itkr- 
We define k\\ as the component of the momentum parallel to 
the local Fermi surface normal, and k± as the component 
orthogonal to the surface normal. 



without affecting the low-energy critical behavior of our 
theory. 

The energy dispersion relative to the true Fermi energy 
at patch a is assumed to be of the form 



?fc =efc°+fc-efc==t'Ffc|| 



—^ = avphn 
2rn 



2m 



, (2.^ 



where vp is the Fermi velocity, k\\ is the component of k 
parallel to the local normal to the Fermi surface, and k± 
is perpendicular to the local Fermi surface normal. In 
the last equality of Eq. (2.8) we have used the fact that 



the Fermi velocity has opposite sign at the two patches. 
We normalize the bosonic field 0^ such that the bare 

fermion-boson interaction vertex Fg ^ ''' [Ki] K2] K^) 
is unity for the model describing the nematic quantum 
phase transition (which we discuss in detail below), and 
assumes the values a = ±1 for the gauge field model, i.e., 

T.^°-0°0/7^ 7^ f> \ _-ra _ \ ^ (nematic model), 
i (Xi; K2; K3) - i - I ^ (gauge model). 

(2.9) 
Finally, the coefficient of the quadratic term in the 
bosonic part of the action is assumed to be of the form 



4-'— »(S) 



Zb — l 



(2.10) 



where po a-nd uq are dimensionful constants with units 
of mass. (Recall that in two dimensions, the density of 
states also has units of mass.) We note that, for a non- 
spherical Fermi surface, mvp does not necessarily equal 
the Fermi momentum. For simplicity, in this work, we 
do not keep track of the renormalization of po by fluctua- 
tions, so that we may set po = to describe the quantum 
critical point. We assume that the bosonic dynamic ex- 
ponent Zf, is in the range 2 < Zf, < 3. To construct a 



sensible limit of large N , the constants po and vq should 
be proportional to TV. However, according to Mross et 
aZ.ES the limit of large- A^ can be safely taken only when 
Zfc— 2 is sent to zero simultaneously, such that the product 
N{zii — 2) remains finite. 

Although the FRG, unlike the field-theoretical renor- 
malization group, does not rely on the presence of a 
small expansion parameter (which indeed is not present 
in the considered problem), it is convenient to express the 
above action in terms of rescaled dimensionless momenta, 
frequencies and fields, to carry out the renormalization 
group procedure. A general discussion of proper scaling 
in mixed Fermi-Bose systems can be found in Ref. l451 
Given an arbitrary momentum scale A, we define dimen- 
sionless fermionic labels Q = {i^,q\\,qi_) by setting 
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(2.11a) 
(2.11b) 
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The corresponding bosonic labels Q = {ie, gy , q±) are de- 
fined in precisely the same way: 



2mvp 
kj_ = Aq_L, 



(2.12a) 
(2.12b) 
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Introducing the rescaled dimensionless fields 
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the Euclidean action of our model can be written as 
SS, V'] = - / E(*^Ae - C^)^§>3., (2.15) 



So[4>] 

"S*! [■0, V', 0] 

where 



{ri, + c^\qi_r'^)(p_Q<t>Q, (2.16) 
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(2.18a) 
(2.18b) 
(2.18c) 
(2.18d) 



and the mixed fermion-boson vertex is the same as before, 







(Qi;Q2;Q3) = r^ '^ ''iK,;Kr,K,) ^r^. 



(2.19) 

If we use the expression i/q = Nra/{2n) for the den- 
sity of states of free fermions in two dimensions, we have 
Co = -/V/(47r). Consequently, setting ta — >■ to describe 
the quantum critical point, our model does not depend 
on any free parameters. 

Our FRG procedure will generate also higher-order 
purely bosonic contributions to the effective action, 
which describe interactions between the boson fields, me- 
diated by the fermions. In an expansion in powers of the 
fields, the lowest-order interaction process is cubic in the 
bosonic fields, 



^3^ = 



^Ki+K2 + K3,Q 
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,Q2,Q3)(I>Q^(I>Q2<PQ3, (2-20) 



with 



r"*'*'*(Qi 



?2,Q3) ^Vp 
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T'f"f"f'(K,,K2,K3). 
(2.21) 



Although for Zf, < 3 this vertex seems to be irrelevant 
by power counting, it turns out that it has a singular 
dependence on the external momenta and frequencies 
and therefore cannot be neglected. Because, within our 
bosonic momentum-transfer cutoff scheme, all vertices 
involving only bosonic external legs a re finite at the ini- 
tial renormalization group scale,'2SIini jt is crucial to keep 
track of the FRG fiow of the vertex T'>">">'{Ki,K2,K:i). 
In this work, we do this by means of a skeleton equation 
relating T'^'^'^{Ki^K2tK^) to the symmetrized fermion 
loop with three external bosonic legs and renormalized 
fermionic propagators, as discussed in Appendix A. An 
explicit evaluation of this vertex is given in Appendix B. 

If we identify A with the renormalization group flow 
parameter which is reduced under the renormalization 
group procedure, the canonical dimensions of all quan- 
tities explicitly appear in the FRG flow equations with 
the above rescaling. To compare the FRG results with 
perturbation theory, it is more convenient, however, not 
to include the canonical dimensions into the definition of 
the vertices. Therefore, we simply choose A = 2mvF in 
the above expressions, so that Ca ~^ 1 and ^a — >• tq = 
po/{2m). Renaming again Q ^ K, our bare action is 



then the sum of the following three terms, 

SS. i^]=- I Y^^lU^- ^kWKarKa, (2-22) 
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where now 



2m' 2m 



The corresponding Gaussian propagators are 

^°^^^ " zw-e^ " tij - ak\\ - kl ' 



F^{K) = h = 
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J^o + colfc^h^-i 
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(2.26) 
(2.27) 



This dimensionless parametrization of our model is what 
is used in the following sections. 



III. FRG FLOW EQUATIONS 

The s tartin g point of our calculation is the Wetterich 
equatiorpSES for the coupled Fermi-Bose model defined 
above, which is an exact FRG flow equation for the gen- 
erating functional V^\%lj,ip,(j)\ of the one-line irreducible 
vertices of our theory. This flow equation describes 
the exact evolution of rA['0i V'l </*] a-s some (now dimen- 
sionless) cutoff parameter A is reduced. By expanding 
Fa ["07 "0,0] in powers of the flelds, we obtain an infinite 
hierarchy of coupled integro-differential equations for the 
one-line irreducible vertices of our model. This hierarchy 
is formally exact, but, in practice, further approximations 
are usually necessary in order to obtain explicit results 
for the vertex functions (see Refs. |36laiid|37|for recent re- 
views). Moreover, the proper choice of the cutoff scheme 
is also very important. 

For our effective low-energy model discussed above, it 
is, in principle, possible to introduce cutoffs in both the 
bosonic and the fermionic sectors and regularize the in- 
verse Gaussian propagators as follows: 



[FoAK)]-' = /; 
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+ Ra{K), 

^1-Ra{k). 



(3.1) 

(3.2) 



For the calculations in the present problem, we find 
it more convenient to introduce a sharp momentum- 
transfer cutoff only in the bosonic sector. Using a sim- 
ilar cutoff procedure, two of us were able to derive the 
exact scaling behavior of the Tomonaga-Luttinger model 
within an FRG approach.!^ We therefore set R\{K) = 
in the fermionic sector, and choose 



Ra{K) 



f^'[e~\\~k±\-A)-i] 



(3.3) 



for the boson cutoff. This leads to the cutoff-dependent 
bare Gaussian propagator 



FoMK)^&{\~k±\-A)h 



fcj 



(3.4) 



which vanishes for |fcj^| < A and equals /^ for |A:j^| > A. 
As there is no cutoff function in the ferniionic sector, all 
purely ferniionic loops already have non- vanishing initial 
values at the beginning of the flow. We see explicitly be- 
low that these are highly singular and need to be treated 
with care. 

The exact hierarchy of FRG flow equations for the one- 
line irreducible vertices of our model can be obtained 
as a special case of the general hierarchy of FRG flow 
equations for mixed Bose-Fermi theories written down 
in Refs. [36] and |40l For our purpose, it is sufficient to 
consider a truncation of this hierarchy which generates, 
after iteration (apart from many other diagrams), the im- 
portant three-loop diagrams identified by Metlitski and 
Sachdev,'221 shown in Fig. [l] Our truncation is character- 



ized by the following three points. 

• On the right-hand side of the flow equations for the 
fermionic and bosonic self-energies, retain only con- 
tributions involving irreducible vertices with three 
external legs. 

• Renormalize all three-legged vertices by triangular 
diagrams involving all combinations of three-legged 
vertices. 

• On the right-hand side of the flow equations for all 
three-legged vertices, approximate the vertex with 
one bosonic and two fermionic external legs by its 
bare value. 

Let us now explicitly give the corresponding FRG flow 
equations. The fermionic self-energy T,°'{K) and bosonic 
self-energy n(^) satisfy the flow equations 
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which are shown graphically in Fig. [3] for a general cutoff 
scheme. Here the scale-dependent bosonic and fermionic 
propagators are 



FiK) 



G°'{K) 
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[i^o,A(if)]-i + n(x)' 
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(3.7) 
(3.8) 



[G^.aW]-'-S"(/0' 

while the corresponding single-scale propagators are 

F{k) = -F\K)dA[Fo,K{K)]-\ (3.9) 

G^{K) = -[G^{K)YdA[GlA{K)]-\ (3.10) 

Note that in the momentum-transfer cutoff scheme 
G" {K) = such that we should omit all diagrams involv- 
ing fermionic single-scale propagators. With a sharp cut- 
off in the bosonic transverse momentum, the full bosonic 



propagator is 



Fa{K) = 



&{\~k±\-A) 



ro + colA^ih"-! + e(|fci| - A)IIa{K) 



, (3.11) 



while the corresponding single-scale propagator is given 
by 



FAiK) 



Si\~k±\-A) 



(3.12) 



ro + coA-"-'' +IU{K) 
The right-hand sides of Eqs. (|3.5|) and (3.6|) also de- 



pend on the one-line irreducible three-point vertex with 
two fermionic and one bosonic external legs, F'^ '^ '^{K+ 
K; K] K), and on the one- line irreducible three-point ver- 
tex with three bosonic external legs, T'^'^'^{K, K' , -K - 
K'), where the superscripts indicate the fields associated 
with the energy-momentum labels. Within our trunca- 
tion, the flow of F''' ''' "^ is determined by the following 
flow equation. 



= ^o^ + ^o- 



(b) 






FIG. 3. (Color online) (a) Graphical representation of the FRG flow equation (3.51 for the fermionic self-energy E" (K) , which is 
represented by a shaded rectangle with one incoming arrow associated with V'ko-j Eind one outgoing arrow associated with tpKcr- 
The shaded triangles represent the three-legged vertex F^ ^ '^{K + K; K; K) with two fermionic and one bosonic external legs. 
The boson propagator is represented by a wavy line. The black dot above the rectangle denotes a derivative with respect to 
the flow parameter, while the slashes attached to the propagators on the right-hand side denote the corresponding single-scale 
propagators, (b) Graphical representation of the FRG flow equation (3.6 1 for the bosonic self-energy Yl{K). The shaded circles 



on the right-hand side represent the symmetrized bosonic three-point vertex. Note that, in the momentum-transfer cutoff 
scheme, all diagrams with a slash on internal fermionic propagators should be omitted. 



OaT 



.0".0" 



\k + k-k-k) 



K' 



K' 



F(k')G°'{K -f K')G°'{K + K + K')+ F{K')G"{K + K')G°'{K + K + K') 

+F{K')G°'{K + K')G°'{K + i? -f K')] r'^°'/^°^(ii' + K'; K; K') 
xr'^°'A°'^(ii: + K;K + K + K'; _if')r'^°'/'°0(^ + K + K';K + K'; K) 
F{k')F{K + K')G°'{K + K + K') + F{K')F{K + k')G"{K + K + K') 

+F{K')F{K + k')G°'{K + K + K')] T'I"'"I'{~K - K' , K, K') 

xr'^>°^(ii' + k;K + K + K'; _i^')r'^°V'°0(^ + k + k';K;k + K'). (3.13) 

A graphical representation of this flow equation is shown in Fig. [4] In the momentum-transfer cutoff scheme, we should 
omit, again, all terms involving the fermionic single-scale propagator. Finally, the flow equation for the symmetrized 
bosonic three-legged vertex is 

= [ [F{K)F{K - Ki)F{k + K2) + F{K)F{K - Ki)F{K -f K2) + F{K)F{K - Ki)F{K + i^a)] 
Jk 

{Ki,K - Ku-k)T'f"f'^{K2,-k - K2, k)T'i"i"i'{-Ki ^k2,~k + ki,k + k2) 



+ f V| [G"(if)G"(i^ + i?i)G"(i^ + ki + k2) -t- G'^iiqc^iK + ki)G''{K + ki + k2) 

+G''{K)G''iK + Xi)G"(if + ki+ k2)]T^'''^"^{K + ki;K; ki)T^'^'^"^{K + ki+k2]K + i^i; X2) 
xr^°V'°0(;^. K + ki+ k2; -ki - X2) + {ki ^ i^a)}, (3.14) 

which is shown graphically in Fig. [5J 



If we ignore the vertex with three bosonic external legs, 
the above system of FRG flow equations has already 
been written down in Ref. I40l (see also Ref. [36| . Al- 
though the above flow equations involve only one-loop 
integrations, the iterative solution of these equations gen- 
erates higher-loop diagrams. In particular, the singular 
three-loop diagrams shown in Fig. [l] are generated as fol- 
lows: The Aslamazov-Larkin type of contribution to the 
bosonic self-energy shown in Fig. [l](b) is generated by 



the last six diagrams on the right-hand side of the flow 
equation for r'*'*'^(Ki, ^2, -^3) shown in Fig. [s] after in- 
tegrating this flow equation over the flow parameter and 
substituting the result into the right-hand side of the 
flow equation for Il{k) shown in Figj3l(b). The singular 
fermionic three- loop diagram in Fig.[lJ(a) is generated by 
substituting the same contribution to r'^'^'^(^i, ^"2, ^3) 
into the flow equation for r'^°'^°''''{K;K';k) shown in 
Fig. |4] After integrating the resulting flow equation 




FIG. 4. (Color online) Graphical representation of the FRG flow equation ([3A3| for the three-legged vertex T''' ''' '^(A'; K'; K) 
with two fermionic and one bosonic external legs. 






+ 




FIG. 5. (Color online) Graphical representation of the FRG flow equation (3.141 for the three-legged boson vertex 



again over A and substituting the resulting vertex cor- 
rection into the right-hand side of the flow equation for 
Y,°'{K) shown in Fig. [31(a), we generate the Methtski- 
Sachdev diagrams shown in Fig. Il](a). In fact, there 
are even further contributions, e.g. using the renor- 
mahzed vertex y^ ^ '^ [K; K'; K) containing the vertex 
T'l>'l>'l>{Ki,K2,K3) as depicted in the second hne of Fig. ^ 
also gives rise to an Aslamazov-Larkin contribution when 
substituted into the diagram on the right-hand side of 
Fig.|3](b). 

Although Eqs. (3.5)-(3.14) form a closed system of 



integro-differential equations for the two-point and three- 
point functions of our model, a direct numerical solu- 
tion of these equations seems to be prohibitively diffi- 
cult, so that further approximations are necessary. As 
a first simplification, we use, below, truncated skele- 
ton e quati ons instead of the FRG flow equations (3.6) 
and (3.141 to determine the bosonic self-energy Il{K) 
and three-point vertex T'^'^'^{Ki, K2, K^). As discussed 



in Appendix A, the Dyson-Schwinger equations of mo- 
tion imply exact skeleton equations relating 11 (^) and 
r'^'^'^{Ki,K2,K3) to the fermionic propagators, the 
three-point vertex F''' ''' "^{K; K'; K) with two fermionic 
and one bosonic external legs and the mixe d fou r-po int 
vertex r'l"f"l"f'{K;K';Ki;K2) (see Eqs. [a^ and (|A3|). 



IV. 



TRUNCATION WITHOUT VERTEX 
CORRECTIONS 



In this section, we show how the known one-loop re- 
sults for the self-energies can be obtained within our FRG 
approach if we neglect vertex corrections. Keeping in 
mind that, in the momentum-transfer cutoff scheme, we 
do not introduce any cutoff in the fermionic sector, the 



(a) 



W 



-H- 





FIG. 6. (Color online) (a) Graphical representation of the 
FRG flow equation (4.5 1 for the fermionic self-energy in the 



momentum-transfer cutoff scheme in the simplest approxima- 
tion where the three-legged vertex is approximated by its bare 
value (represented by a black dot), (b) Truncated skeleton 
equation (4.41 for the bosonic self-energy. 



scale-dependent fermionic propagator is 



GliK) 



tuj-^^-[j:i{K)-j:i{o)Y 



(4.1) 



where S^(0) is the self-energy at the renormaUzed flow- 
ing Fermi surface. The subtraction of S^(0) is necessary 
because, by assumption, we have expanded the wave vec- 
tor at the renormalized Fermi surface of the underlying 
model. Given the cutoff-dependent self-energy 'Sa{K), 
we may define the Fermi surface for a given value of the 
cutoff parameter A via 



efep +^A{iO,kp) = fi, 



(4.2) 



which, for A — > 0, reduces to the definition of the true 
Fermi surface. Hence, we may write 

efe=+fc - M = efcj+fc - Efcp - SA(iO, fc?) = Cfc - Sa(0). 

(4.3) 
Following Refs. l43l l44l and |47l we now use the skele- 



ton equation ( Al I to determine the flowing bosonic self- 
energy. Using the fact that (Fq)^ = 1, the scale- 



dependent bosonic self-energy is thus given by 

nA(^) = ^ / E ^A(i^) [GliK + K) + Gl{K - K)] , 

(4.4) 
while the fermionic self-energy satisfies the FRG flow 
equation 



dAj:i{K) = / Fa{k)gi{k + K). 



(4.5) 



K 



A graphical representation of Eqs. (4.4) and (4.5 1 is 



shown in Fi g.[6l Actually, the integral on the right-hand 
side of Eq. (4.5 1 is ultraviolet divergent for our model, 



but since wc do not keep track of the shape of the renor- 
malized Fermi surface, we can instead consider the sub- 
tracted self-energy, 



/^l{K) = Y.l{K)-Y.l{Q), _ ^ (4.6) 
which appears in our cutoff-dependent fermionic propa- 
gator (4.1). The subtracted self-energy then satisfies the 



FRG flow equation 



df^^HK) = f hm [G%iK + K)- GliK)] 
Jk 



(4.7) 



To obtain an approximate solution of this flow equa- 
tion, we expand the subtracted self-energy for small mo- 
menta and frequencies. 



^liK) 



-izi' 



1)^0. + (Z^i - 1)^1 + 



(4.8) 



where we have used the fact that the fe-dcpendence of the 
self-energy appears only in the combination i^^ . Because 
the cutoff. A, regularizes the singularity of the bare in- 
teraction, the self-energy is analytic for small momenta 
and frequencies and can hence be expanded into a Taylor 
series. The corresponding low-energy approximation for 
the fermion propagator is 



GliK) 



ZlHu - Zl^il luj - (Za/Za)C^ 



(4.9) 



Substituting Eq. (4.9) into Eq. (|4.4[), it is convenient to 
first perform the integration over fcii ,^^ which can be done 
using the residue theorem. 



nA(i^) = ^ZAZA/%E 



27r 



^ iCj- (ZA/^A)(afc|| +kl + 2k^k^) 



duj 
2^ 



eiaioj + oj)) - Qiauj)] + iK ■ 



-K). 
(4.10) 



The integration over u is now trivial. 



— [e(a(a. + a;))-e(aw)] = — , 



(4.11) 
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so that 



n.m-lz..zJ'-t^ 



iCj - {Zf,/Zf,)(ak\\ +k\+ 2kj_kj_) 



(K -^ -K) 



(4.12) 



r 



Note that the fcj_ integral is stiU ultraviolet divergent. 
Following Mross et a/.p^lwe regularize the divergence by 
symmetrizing the integrand with respect to A;j_ o — A:_l, 
so that the expression in the square braces vanishes as 
l/k\ for large k^. The fcj^ integration can thus again be 
done using the method of residues, with the result 



nA(w,fc^) =6a 



fc±| 



where 



h ^72 



(4.13) 



(4.14) 



Note that this is independent of fcii , the mathematical 
reason being that the term ak\\ in the denominator of 
Eq. (4.12 1 can be eliminated by means of a simple shift of 
the integration variable k±^ . We show shortly that at one- 
loop level, the self-energy S^ [ K) is actually independent 
of fc, so that Za = 1 and Eq. (4.131 reduces to 



where 



nA(w,fcj 



&0 



= &0 

TV 
4^' 



w 



\ki.y 



(4.15) 



(4.16) 



in agreement with Metlitski and Sachdev.--^- It should be 
noted that the bosonic self-energy does not renormalize 
the exponent Zb which characterizes the momentum de- 
pendence of the bare boson p ropag ator.!^ 

Next, we substitute Eq. (4.13) into our expression 



(3.12) for the bosonic single-scale propagator and obtain 



Fa{K) 



5{\k^\-K) 



»'o 



(4.17) 



coA^i.-i-hfoAl'I'l/A' 

Substituting this expression into the FRG flow equation 
(4.7) for the subtracted self-energy, we may perforin all 



integrations on the right-hand side and obtain 



ttoa 27r 



1 



^aI^I/a 

To -I- Co A^"-! 



(4.18) 
Relating both Za and Za to flowing anomalous dimen- 



sions via 



i-jK = A9a In Ziy, fJA = Ad A In Za , 



(4.19) 



the corresponding flowing "frequency" anomalous dimen- 
sion 77a is given by"^^ 

d 

r]A = AZa lim —--dAA%{iuj, fc = 0) 



>o d{iuj) 
AZaZa 



27r2(ro -l-coA^f'-i)' 



(4.20) 



Keeping in mind that the right-hand side of the flow 
equation (4.18) for the self-energy depends only on uj 
and not on fe, we conclude that Za = 1 in this approxi- 
mation, so that bA — bo. At the quantum critical point, 
tq = 0, our expression (4.20) for the flowing frequency 



anomalous dimension hence reduces to 

Za 

"^^ 27r2coA^'.-2 ' 

so that Za satisfies the flow equation 



AOaZa = vaZa 



72 



(4.21) 



(4.22) 



27r2coA^^-2' 
This differential equation can be easily integrated to give 

Zao 



Za 



1 



^Aq 






2-Zi, 



Al^') 



(4.23) 



Assuming zt > 2, we see that, for A ^ 0, the wave- 
function rcnormalization factor vanishes as 



ZA^27r^co{zb-2)A''-\ 



(4.24) 



implying that, at the quantum critical point, the fermion 
field has the frequency anomalous dimension 



V 



lim 

A-s-O 



Ad In Za 
dA 



Zb 



(4.25) 



Because, for Za = 1, the right-hand side of Eq. (4.18) 



depends on the flow parameter A only via the explicit A- 
dependence shown, we can obtain the flowing subtracted 
self-energy A'j!^{K) for all frequencies by simply integrat- 
ing both sides of Eq. (4.181 over A, 



A%iK) - Al^{K) = - dXdxAtiK) 

J A 



2n^bo 



1 



boH/X 



To + CqX'"'- 



(4.26) 



For Zb > 2, the integral is ultraviolet convergent, so that 
we may take the limit Aq — )■ c« where A^ (K) — >■ 0. At 
the quantum critical point, tq = 0, the dependence of the 
integral on uj can be scaled out and we obtain for A — >■ 



jimAKK) 

A— i-O 



-isgn a; — 
Oo 



5oM 
Co 



2/zi 



where 



Oo 



2tt^ Zb 



"•^ 2.1+2/z, 



47r sin(27r/zb) 



(4.27) 



(4.28) 



in agreement with previous workPSESIMSIl ^]^q corre- 
sponding one-loop corrected fermion propagator is thus 



G°{K) 



1 



1 



bo\uj\ 

Co 



2/zb 



2/26 



i^i^)Z{t^) '-<^h-kl 



,(4.29) 



where we have used the fact that for Zb > 2 the Mat- 
subara frequency ito is smaU compared with the contri- 
bution from the self-energy at low energies. Comparing 
Eq. (4.29) with the general scaling form (1.2), we con- 
clude that z = Zb/2 and r/^ = within the one-loop 
approximation. We note that this is in agreement with 
the one-loop results of previous works based on the field 
theoretical rcnormalization group . ^^ ^^ 29 32 



V. VERTEX CORRECTIONS 

In this section, we use the hierarchy of FRG flow equa- 
tions given in Sec. |III| to estimate the effect of vertex cor- 
rections on the low-frequency behavior of the fermionic 
self-energy, S]^(iir). We do not attempt to calculate the 
entire momentum and frequency dependence of E^(i4r), 
but focus on the flow of the two renormalization con- 
stants Z\ and Z\ deflned via the low-energy expansion 



(4.8) of the flowing self-energy, and on the associated 



anomalous dimensions 77A and ?7a defined via the loga- 
rithmic derivatives of Z\ and Za with respect to the flow 



parameter (see Eq. (4.19)). Note that the deflnition (4.9) 
implies 



1 — lim 



Z7^ = 1 + lim 



d{iuj) 

9EX(ia; = 0,fc) 



d^i 



(5.1) 
(5.2) 
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while the definition (4.191 of the flowing anomalous di- 



mensions, ?7a and tja, allows us to relate these quantities 
directly to the derivative of the self-energy with respect 
to the flow parameter 



?7a = A9a lu Za = ZaA lim 



d 



c^^o d{i(jj) 



dAEl{lLO,k^0), 



(5.3) 



fJA = A^A In Za = -^aA hm --dA^%{tuj = 0, fc) 



= — ZaA lim 



d 



fc-i-o 9(afcj| 



-dA^l{lOJ = 0,k), 



(5.4) 



where, in the last line, we have used the fact that the 
momentum dependence of the self-energy appears only 
in the combination ^^ = ak^^ + k'^. In the one-loop ap- 
proximation, liniA-i-o Va = 0, but in this section we show 
that vertex corrections lead to a finite value of this limit. 
In Sec. VD we further show that liuiA-i-o ^A can be iden- 
tified with the anomalous dimension 77^ of the fermion 
field defined via Eq. ( 1.2 1; moreover we show how to ex- 
press the fermionic dynamic exponent z in terms of rjA 
and 77a- 



Truncated flow equations and skeleton 
equations 



Using the momentum-transfer cutoff scheme in combi- 
nation with the truncation strategy discussed in the third 



paragraph of Sec. Ill we obtain the fermionic self-energy 
from 



where the three-point vertex with one bosonic and two fermionic legs is determined by 



(9AE"(i^)= / F{K)G'^{K + K)Tf 
Jk 



(5.5) 



gjj^i'''i'°'4'(^K + K;K;K) = iT^f / F{K')G°'{K + K')G°'{K + K + K') 



K' 



{T^f [ F{K')F{K + K') [G^iK + K + K') + G^'iK - K')] r'*"'"''{K, K' , -K - K'). 

JK' 



(5.6) 



r 



Eq. (5.51 can be obtained from the more general flow three-legged vertex in the second line of Eq. (5.6) is ap- 



equation (3.5) by simply omitting the contribution in- 



volving the fermionic single-scale propagator, while in 
Eq. (5.6) we have, in addition, replaced the flowing ver- 



tices T'^ ^'''P{K + K; K; K) on the right-hand side of the 
flow equation by their bare values. The purely bosonic 
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(a) 



(6) 



= ^^3_ 



(c) 



y^ 




FIG. 7. (Color online) (a) Graphical representation of the 
FRG flow equation (5.5 1 for the fermionic self-energy in the 
momentum-transfer cutoff scheme. This flow equation can 



be obtained from the more general flow equation (3.5 1 by 



omitting the contribution involving the fermionic single-scale 
propagator, (b) Truncated flow equation for the three-legged 
vertex r'''" '''''''' {K + K; K; K) in the momentum-transfer cut- 
off scheme, which can be obtained from the more general flow 
equation (3.131 shown in Fig.lllby neglecting the contribution 
from the fermionic single-scale propagator and approximating 
the vertices T'''"' 'I'"' "t- [K + K; K; K) on the right-hand side by 
the bare vertices, (c) Graphical representation of the approxi- 
mate expression (5.71 for the bosonic three-legged vertex used 
in our calculation. Note that this expression can be obtained 
from the exact skeleton equation ( A3 1 derived in Appendix A 



by neglecting the irreducible four-point vertex and then mak- 
ing the same approximations as in (b). 



proximated by 

G°'{K)G"{K + Ki)G°'{K + /?i + K2) 



K 



G"{K)G°'{K + Ki)G°'{K - K2) 



K 



, (5.7) 



which is derived from the exact skeleton equation by 
making the same approximations as in the derivation 
of Eqs. (5.5) and (5.6). A graphical representation of 
Eqs. (5.5)-(5.7) is shown in Fig. It] To obtain the cor- 
rections to the one-loop approximation, it is sufficient to 
approximate all propagators in Eqs. (5.5)~(5.7) by the 



one-loop results given in Sec. IV i.e., we approximate on 



the right-hand sides 
Zk 



Gl{K) = 



Fk{K) = 



e(lfc^l-A) 

6{\k^\-K) 



(5.8) 
,(5.9) 
(5.10) 



Z\ is the one-loop result for the flowing wave- 



Here, 

function renormalization factor given in Eq. (4.23) and 



we have used the fact that Zf^ — 1 within the one-loop ap- 
proximation. We thus arrive at a system of flow equations 
for the momentum- and frequency-dependent two-point 
and three-point functions. Having determined the right- 
hand side of the flow equation (5.5) for the self-energy. 



we can substitute the result into Eqs. (5.3) and (5.4) and 



obtain for the flowing anomalous dimension related to 
the frequency dependence of the self-energy 

,7A = ^aA lim -^ / F^{K)[Gl{K + K) - Gl{K)] 



x[r 



A 



+ 2ZaA 






K 



F^{K)Gl{K)Tl ^ ^{K-Q-K) 



-^Ip^ Ip^ (f) f 



dTl^'^(K + K:K;K) 

X lim jTT — ^ — '■ . 

K^Q d{iuj) 



(5.11) 



and for the corresponding anomalous dimension associ- 
ated with the momentum dependence of the self-energy, 

77A = -ZaA lim ^J^^ f F^{K)[Gl{K + K) - Gl{K)] 

IK 



'k^o d{ak\\ 



- 2ZaA 



x[rf'^°*(K;0;K)]2 
K{K)Gl{K)Tf^''^{K-Q-K) 

K 

drf^°^{K + K;K;K) 
a:->o d[ak\\) 



(5.12) 



In the first terms on the right-hand sides of these expres- 
sions we have used the same regularization as in Eq. (4.7 1 



B. Bosonic three-legged vertex 



In order to evaluate the self-energy from Eq. (5.5 1, 
we should calculate the mixed fermion-boson vertex 
P^°^°0^^ _^ K;K;K) by integrating the flow equation 
(5.6), which in turn depends on the purely bosonic three- 

K2). Fortunately, the 



legged vertex T'l"l"l'{Ki,K2,-Ki 

integrations in the truncated skeleton equation (5.7) can 



be carried out exactly in our model, as shown in Ap- 
pendix B. The result can be written as 

r***(K,,K2,-K,-K2) 

= 2\NY,{nfL'^{K^,k2,-k^ - K2), (5.13) 

a 

where the symmetrized fermion loop with three external 
bosonic legs and one-loop renormalized fermionic propa- 
gators is given by 
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L^iK,,K2,-Ki^K2) = 



1 



1 



An \kj^i 



1 

k±2 



-10 te)+-^© (&)-(' 



wi + u;2 



UJ2 






t(&-s) -(^• 



- fc±2)2 



(5.14) 



Note that this function represents a rather complicated 
momentum- and frequency-dependent effective interac- 
tion between the bosonic fluctuations, mediated by the 
fermions. Obviously, this function cannot simply be ap- 
proximated by a constant, which is assumed to be pos- 
sible in th e Her tz-Millis approach to quantum critical 
phenomenaliSEil 



C. Three-legged boson-fermion vertex 

We have now calculated all functions appearing on the 
right-hand side of the FRG flow equation (5.6) for the 



three-legged vertex with one bosonic and two fermionic 
external legs, so that we may next integrate this equation 
over the flow parameter A. Let us begin by evaluating 
the first term on the right-hand side of Eq. (5.6), 



dj^r'''''''''''l'{K + K;K;K) 



(1) 



(r? 



Q\3 



F{K') 



K' 



cG'^iK + K')G"{K + K + K'), 



(5.15) 



corresponding to the first diagram on the right-hand side 
of Fig. [71(b). The integrations in Eq. (5.15) can be ex- 
plicitly carried out, with the result 



KdAr'^'^''''^'l'(K + K;K;K) 



(1) - 



i^^Y Za 



A2 



(27r)2 6o 
i sgn(ci; + oj) In ( 1 + 



iu) - ZA_[ak\\ +k\ + 2k±{k± + A) 



roA + coA^ 



i sgn(aj) In 1 + 



+ (A ^ -A) 
bo\uj 



roA + coA^f 



(5.16) 



For the evaluation of the flowing anomalous dimen- 
sions in Eqs. (5.11) and (5.12), we need only the 
vertex F''' '^ '^{K;0;K) at vanishing fermionic energy- 
momentum, as well as the derivatives of F''' '^ '^{K + 
K; K; K) with respect to the components oi K dX K = Q. 
From Eq. ( 5.16[ ), we see that the contribution of the first 
diagram in Fig. [71(b) to the flow of these quantities can 
be written as 



^^^pV;°V'°0(^.O;if)(l) = 



A^A 



5pV'°'A°0(X;if;ii')(l) 



d{iuj) 



1 



r(l") 



A^A 






d{akii 



^A 


/ Cj /c|| k^\ 
\ A^f' ' A2 ' A ) 


,(5.17) 


K=Q 








h ~kA 
A2 ' Ay' 


(5.18) 




= 0, 


(5.19) 



K=0 



where we have defined the dimensionless scaling functions 



iA\<^,q\\,q±) 



7A . .-X, /, , bale] 

--—-I sgn e In 1 + —— 

(27r)2feo V rA + Co 



1 



ie-jA[aq\\ + qj_ + 2qj_] 



(qi 



-q^) 



, (5.20) 



and 



r(l") 



^A (e.^lh^-LJ 



7A 



(27r)2(rA + co)2l + ^|e-| 



ie-7A[a<?|i + q'i + 2gj_] 



(g_L -^ ~q±) 



(5.21) 



Here, ta — tqA^^^'' vanishes at the quantum critical 
point, and wc have introduced the dimensionless param- 
eter 



7A 



Za 



27r2co(2;fc-2) 



A^f'-2 i-A^i>-2r 



_2r27r^Co(zt,-2) 



ZAn 



-A^-] 



, (5.22) 



which approaches the limit 2TT^co{zb — 2) ex N{zb — 2) for 
A^O. 

Next, consider the contribution of the last two dia- 
grams on the right-hand side of Fig. [7|(b) to the flow of 
the three-legged boson-fermion vertex, corresponding to 
the second line in Eq. (5.6), 
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^^P^°^°,/,(|^ + K; K; K)'-^'^ = -{Tof f F{K')F{K + K') [G"(if + K + K') + G"{K - K') 

JK' 



:T't"t"t'{K,K',-K-K') 



(5.23) 



Due to the rather complicated form o f the vertex T^^'^jK, K', —K — K') given in Eqs. (5.131 and (5.141, the evaluation 



of the right-hand side of Eq. (5.23) is quite involved. The k',, integration can still be performed by means of the 



residue theorem, while the k'j_ integration is trivial due to the i5-function in the single-scale propagator. After these 
integrations, we obtain 



g^pV'°V'°0(^ -K K; K- K)^^^ = -iaNJS? 



Q\2 foa 



(27r)2 



du' 



1 



27r ro + coA^"-! +ho\uj'\/K 



^F{Q + Q' ,k^ + sK)J{oj,k^,oj' ,sK) Y^ {T^'f 



s=± 



a'=± 



e(Im(zi)) - e(Im(z3)) e(Im(z2)) - e(Im(z3)) 



{zi - zs)^ - zl 



{Z2 - Zs)^ - zl 



(5.24) 



where 



J(a;,fcj_,a;',fc^) 



1 1 



^^'Sj+'^'^. 



{u} + uj')Q 



kl 



(5.25) 



and 



Zl 



-{k\\+k{)-a[K + s{k^ + k^)y 

id)' + id) + ILU 



Zk 



Z2 = fc|| + ck(A — sfc_L) + a- 

A 



zz 

Zi 



' fc|| a' 

skj_ Zf^ 
A(A + sfcj_). 



sk± 



Za 

id) 



(5.26a) 
(5.26b) 

(5.26c) 
(5.26d) 



For the calculation of t^a and t^a in Eqs. ( |5.11 ) and (5.12) 
we, again, need only the vertex T^ ^ '^{K; 0; K) and the 
derivatives of r'^°'^°'^(i? + K; K; K) with respect to the 
components of the fermionic label K at K = 0. In anal- 
ogy with Eqs. (5.17)~([5.19 1, we therefore define 



^^^pV;°V°0(^.O;if)(2) 



r(2) 



U) /cy fcj_ 

A^' A2'T 



,(5.27) 



A^A 



arV'°'»(ii:;if;ii:)(2) 



d{iu)) 



1 .(2a.) 



Kds 



5P^°^°0(^.^.^)(2) 

9(afcj|) 



-i.' 



{2k) 



K=0 

dl fc|| fcj_ 

A^' A2'T 



K=0 
OJ fc|| fcj_ 

A^' A2'T 



(5.28) 



(5.29) 



From now on we focus on the Ising-nematic instabil- 
ity and explicitly set Fq = 1. In principle, the functions 

-^A (e,9|h9-i), -^A ""(e, 911,9-1 ), and l';^'^\e,q\\,q±_) can be 



calculated analytically by performing the remaining fre- 
quency integration in Eq. (5.24). However, the result is 



rather cumbersome and not very transparent. To sim- 
plify the calculation, let us assume that the arguments 
e, q\\ , and q± of these functions are all small compared 
with unity. Then the resulting expressions simplify and 
we obtain the approximate expressions 



-^A (e,9i|,9±) 



^A (e>9i|,9±) 



7A|e|^ 



2(2^)3(rA + co)2 



1 



\qj-\ -w 
\'i±\ 



w 

«7Ag|||£| 
e^ + 7a 9y 



r(2fc) 



4 (e,9lh9-L) 



(k-^|-u;)2^^ 

7A^^ 

(27r)3(rA-Hco)2 

e'(e'-7ig^) 

-^sgn(e-)k-J ^.,^^,_,^, 

isgnie)jl\q^\NeHe^-llqj) 
(27r)3(rA + co)2 (£-2+72^-2)2 



(5.30) 



(5.31) 



, (5.32) 



where, in Eq. (5.30), the complex quantity w is defined 
by 



w 



bo 



2(rA + Co) 



(|e| - Ja7Ag|tsgn(e)) 



(5.33) 



Let us now combine the contributions from all three 
diagrams on the right-hand side of Fig. [71(b). To be con- 
sistent with the approximations made in the derivation of 
Eqs. ( 5.30 )-( 5.32), we should also expand the right-hand 



sides of /\ {e,q\\,q±) and I}^'^ {e,qi\,q±_) in Eqs. (5.20) 



and (5.21) for small e. We therefore define 
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^A(e,'7||,g_L) = /a (e,g||,g±) + 1^ (e,g||,g±) 

7A J ie 



le 



(27r)2(rA + Co) I ie - ^A[aq\\ + ?! + 2g_L] «£ - 7A[ag|| + ?! - 2(7j.] 



TV 



47r(r'A + Co) 



\q±\-w {\qj_\~w) 



q±\\e\ i^/l^-Ll 



7A^A J l« 



(27r)2(rA + co)2 I ie - 7A[agj| + 9^ + 29_l] «e - 7A[ag|| + ?! - 2g_L] 



TV 
"27r6A 



a7Ag| 



|e| Z6|e|(62 - 72g2) 



e" + 7a 9' 



19^1 



i^^+lA<in) 



2;72\2 



2;72\ 



-^a (e.9lh9-L) = /a (e.^lh^i) = 



2 =2\2 



(27r)3(rA + co)2 (6-2+72^-2) 



(5.34a) 



(5.34b) 
(5.34c) 



Finally, we integrate over the flow parameter A, and ob- 
tain the following expressions for the three-legged boson- 
fermion vertex, 



-^ip^ ^^ <p 



{K-Q-K)^TA{e,q\\,qi:), (5.35a) 



dT'l''''i''''*'{K;K]K) 



d{iuj) 
dT'!''''t''''>'{k;K;K) 



d{ak\\ 



K=0 



K=0 



^fX(e-g||,gi), (5.35b) 



A2 



fi(e-,gj|,g-_i), (5.35c) 



where 

f^ ds 

^A{e,q\\,q±) = 1+ / — ^A/s (s^'e, s^gy, sg_L) , 
Ja/Ao s 

(5.36a) 

JA/Ao 

(5.36b) 

f ^(e, 911 ,qi_)= f ds sll,^ (s^-'e, s'm , sqA . (5.36c) 



A/Ao 



Recall that in deriving these expressions we have assumed 
that lei < 1. 



D. Fermionic anomalous dimension and dynamic 
exponent 

We are now ready to calculate the anomalous dimen- 
sions 77 — limA-i-o "Ha and f] = liniA^o ^A at th e quan - 
tum critical point. We therefore substitute Eqs. (|5.35aP - 
( |5.35c l into our general relations (5.11 ) and ( 5.12[ ) forThe 
flowing anomalous dimensions 77A and f]A and, introduc- 
ing the integration variables p = ^AO^q\\ and y = e, we 
obtain for the flowing anomalous dimension associated 
with the frequency dependence of the self-energy. 



??A 



7A ,. d 

--— lim —7 — - 

27r w-i-o d{iijj) 

TA [^° dy_ 
' ^ J-x„ 27r 



27r 



-A, 
dp 

-^ 27r 



dp 

2^ 

1 



1 



rA + co + 6o|y| 

1 



T-A + co + fco|y| 



1 



1 



W-P-IA 



iy + iuj - p - jA iy-p-jA 



rA[y,^,i)rt{y,^,i 

ajA J V "7A 



Fa y,^-,l 

a7A 



(5.37) 
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and for the corresponding anomalous dimension associated with the momentum dependence of the self-energy, 



7A ,. d 
2tt fe-i-o ok 



-A 
Ao j^, /"oo 



2n 



° dy 

2tt 

dp 

2ti 



dp 

2ti 



1 



Co 



rK 



1 



%\y\ 



ly 



1 



1 



iy - p - k - jA 



ly-p-^A 



P 

a-lh 



rA ( 2/, T^— ,1 



1 
p 



7A 



rA y 



a7A' 



1 r 



A V 



ajA' 



(5.38) 



Here, Aq is an ultraviolet cutoff of the order of unity which takes into account that our expressions ( 5.34a )-( 5.34c) 



which we use to calculate the vertices in Eqs. ( 5.36a )-( 5.36c I are only valid for small frequencies. 

Consider now the limit A — >■ 0. At the quantum critical point we may then set ta — > 0. For simplicity, we also 
set Co = bo = N/{4:Tt). To make progress analytically, let us further assume that the parameter 7 = limA -»o 7A = 
2Tr'^co{zi, — 2) = ^iV(z(, — 2) is small compared with unity. To leading order in zt, — 2 the s-integrations in Eqs. ( 5.36a)- 



(5.36c I can then be carried out analytically, with the result 



lim Fa 

A-J-O 



y, 



p 

a-jA 



1=1 



7 
(27r)2co 



^y_ 

iy-p 
iyN 



In 



7 



1- 



In 



ly-p- 

27 

2isgn(y) -| 



ln[l 



iy+p 1 

2isgn(K)J 



lim F 

A^O 



A y, 



lim ft 

A-s-O ^^ 



y, 



p 

a7A' 

P 
a-jA 



1 = 



27rco 

7 [ \y\ ^ N 

(27r)2co\ iy-p 2ttco 
7^JV iy\y\{y'^ - p^) 

(27r)3c2 3(2/2 +p2)2 ■ 



2zsgn(?/) - {iy+p) 

\y\p iy\y\{y'^ 



iy+p 



y2 ^ p2 3(y2 + p2) 



P^ 

2^2 



(5.39) 
(5.40) 
(5.41) 



Note that these vertex functions describe the renormal- 
ized effective interaction between two fermions and one 
boson; clearly, this interaction has a rather complicated 
dependence on momenta and frequencies and cannot 
be approximated by a constant. If we now substitute 



Eqs. ( 5.39 1-( 5.41) into Eqs. (5.371 and (5.381, we may 



perform the p-integration analytically by means of the 
method of residues. Note that the first term in the curly 
braces of Eq. ( 5.39 1 and the first term in the curly braces 



of Eq. (5.40) do not contribute to the integrals, because 



we may close the integration contour in a half plane where 
the integrand is analytic. These terms arise from the first 
diagram on the right-hand side of Fig. [71(b), so that the 
three-boson vertex is crucial to obtain the leading cor- 
rections to the anomalous dimensions. We finally arrive 
at the following estimates for the anomalous dimensions 
at the nematic quantum critical point, 

r, = Z6 - 2 + i^^^^aXo) + 0{{z, - 2)3), (5.42) 



V^ ^'\^^^ C{\o)+0{{z,-2f), 
where the cutoff functions are given by 



C{\o) 



C{\o) 



'" dy 



l + y 
'° dy 



l + y 



1 I ^ I y^^y 

6 1 — y (1 — y)2 
1,1, ylny 



(5.43) 



,(5.44) 



_6 l-y {l-y)\ 
= C(Ao)-ln(l + Ao). (5.45) 

A plot of the cutoff functions is shown in Fig. [8J Note 




0.5 



FIG. 8. (Color online) Plot of the cuto ff functions C(Ao) and 
C'(Ao) defined in Eqs. |5.44| and ([s^isi. 



that these functions depend logarithmically on the ultra- 
violet cutoff, Aq. This is due to the fact that we have 
assumed that the rescaled frequencies and momenta in 
Eqs. ( |5.30[ )-( |5.32 | are small compared with unity in the 
evaluation of the vertex-correction diagrams shown in 
Fig. 0(b). By retaining the full momentum- and fre- 
quency dependence on the right-hand side of the flow 
equation ( |5.24 ), one obtains ultraviolet convergent re- 
sults from our ERG approach. However, as the cutoff 
dependence in our final result, Eqs. (5.42 1 and (5.43), is 



only logarithmic, we have not attempted to carry out this 
calculation which requires substantial numerical effort. 
Instead, we make the reasonable cutoff choice Aq = 1 
and find, within the accuracy of our calculation, that 
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C w 1.27 and C w 0.58. 

Given our result for the anomalous dimensions ry and 77, 
we may relate these to the anomalous dimension rj^ of the 
fermion field and the fermionic dynamic exponent, z. Let 
us therefore recall that, at the quantum critical point, the 
retarded single-particle Green function assumes for small 
frequencies and momenta the following scaling form. 



G"(w + iO+,fe)oc 



1 



[^^sgnH|w|iA-e^ 



l-')^/2' 



(5.46) 



where Ai^ is some dimensionful constant with positive 
imaginary part and real part depending on the sign of oj, 
see Eq. (1.2). For u = k± = 0, this implies 



G"(iO+,/c||,0)(xfc^^+''*/^ 



(5.47) 



On the other hand, from the definition of 77A in Eq. (4.19) 
we see that Z\ oc A^ ex fci? , where we have used the fact 
that fell scales as A^. Hence, 



G"(iO+,fc|[,0) cxZA/fc|| ocfc. 



-l+fj/2 



(5.48) 



Comparing this with Eq. (5.47), we conclude that 



Vi> = 1- 



Next, setting A; = in Eq. (5.46), we find 



G"(a; + iO+, 0) oc cJ-(l-'''*/2)/^ 



(5.49) 



(5.50) 



On the other hand, from the definition of rj\ in Eq. (4.19 ) 
we infer that Za oc A'' ex aj''/(^^) (using the fact that lu 
scales as A^^), so that 



G"(0,a; + iO+) oc Ziy/uo oc w-(i-''/(2^)). 



(5.51) 



Comparing this with Eq. (5.50), we see that 1 — r//(2z) 
(1 - 'n^/2)/z, or 



z = 1 



rj — rj 



Writing Eq. (5.42) as 77 = z;, — 2 + (5?7 with 



5^1 - ^^^-^G(Ao) + 0{{z, - 2)3), 



(5.52) 



(5.53) 



we see that Eq. (5.52) can also be written as 



Zb , Srj-f] 

^ Zb {zb - 2f 
2 4 



G(Ao)-G(Ao) +0((zb-2)3). 

(5.54) 



According to the above calculation, the scaling relation 
z = Zfc/2 for the fermionic dynamic exponent acquires 
a correction of order [zb — 2)^, such that the effective 



theory appears to have two different time scales. How- 
ever, according to Metlitski and Sachdev,'221 the strong 
correlations imply that z ~ Zb/2 should be satisfied ex- 
actly. This discrepancy might be due to the fact that 
our truncation of the FRG flow equations introduces ap- 
proximations which violate the general scaling theory. 
In fact, this seems to be true also in the case of the 
loop ex pan sion. In the previous works of Metlitski and 
Sachdev^^l and Mross et a/.,'^the three-loop correction 
to the fermionic anomalous dimension rj^ is calculated 
via the momentum dependence of the self-energy, and 
according to scaling theory it is argued that a corre- 
sponding correction should appear if one calculates rj^ via 
the frequency dependence of the self-energy. However, as 
pointed out in Ref . 1311 the correction in the frequency de- 
pendence of the self-energy (the term (5r/ in our notation) 
appears only beyond the three-loop order. The FRG is 
not based on an expansion in powers of loops, so that dia- 
grams beyond three loops are included in our truncation. 



VI. SUMMARY AND CONCLUSIONS 

In this work, we have used a functional renormalization 
group approach to calculate the anomalous dimension 77^ 
of the fermion field and the fermionic dynamic exponent 
z of an effective low-energy field theory describing the 
Ising-nematic quantum critical point in two-dimensional 
metals. In the limit N{zb — 2) ^ 1 (where N is the 
number of fermion flavors and Zb is the bosonic dynamic 
exponent), we have been able to explicitly calculate the 
fermionic anomalous dimension 77^, of the system, with 
the result 



77V,«0.3(zfc-2)'+O((zb-2)3). 



(6.1) 



If we extrapolate this expression to the physically rele- 
vant case, Zb = 3, we obtain rj^ sa 0.3, which is larger 
than the estimate 770 « 0.068 given by Metlitski and 
Sachdev,^^ but still smaller than the estimate 77^ « 0.6 
obtained from an extrapolation of the corresponding ex- 
pression given by Mross et al.^^ Given the different types 
of approximations, it is not surprising that different val- 
ues for the exponent rj^p are found. This fact shows that 
the recent calculations done so far are not yet fully un- 
der control in the physically relevant case oi Zb = i and 
N = 2. It is however reassuring that our method also 
finds a finite value for 77^. In particular, we note that in 
the previous calculations based on the field-theoretical 
renormalization group, the values of the exponent 77^ are 
obtained within a loop-expansion truncated at the three- 
loop order, while in our FRG approach the truncation 
does not rely on a loop-expansion, so that a certain class 
of diagrams beyond the three-loop level are effectively 
re-summed to higher order within our truncation. On 
the other hand, at the lowest order, namely at the one- 
loop level in the field theoretical RG and neglecting the 
vertex corrections in the FRG, the results obtained us- 
ing the FRG and other renormalization group methods 
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FIG. 9. (Color online) Exact skeleton equation relating the 
bosonic self-energy to the exact fermionic propagators and 
the thre e-leg ged vertex with two fermion and one boson leg, 
see Eq. ( Al I. The black dot denotes the bare vertex with two 



fermionic and one bosonic external leg given in Eq. (2.9 1. 



all coincide with the well-known results obtained w ithin 
the RPA. Finally, in contrast to previous works ,'22E11 
we have explicitly shown that both the frequency and 
the momentum dependence of the self-energy give rise to 
anomalous corrections to the one-loop result. While our 
calculations in principle lead to a small correction term 
to the fermionic dynamic critical exponent, z, the scal- 
ing theory of Metlitski and Sachdev^ implies the exact 
identity z = z-^j^. 

The calculations presented in this work can be ex- 
tended in several directions. Because our FRG approach 
does not rely on the smallness of the parameter A^(zfc — 2), 
with some numerical effort it should be possible to ex- 
tract the anomalous dimension 77^ for general TV and z;,— 2 
from Eqs. (5.371 and (5.38). In this case the vertices ap- 



pearing in these expressions cannot be calculated analyt- 
ically but must be represented as one-dimensional inte- 
grals, so that the evaluation of the anomalous dimensions 
in Eqs. (5.37) and (5.38) requires rather complicated nu- 
merical integrations, which is beyond the scope of this 
work. In principle our FRG flow equations also allow 
us to calculate the entire frequency-dependence of the 
self-energies, but this seems to be numerically even more 
expensive. 

In cont rast t o the strategy adopted in the Hertz-Millis 
approach jlSin] in the present problem it is not possible to 
integrate over the fermionic degrees of freedom to obtain 
an effective bosonic theory with regular vertices. We have 
therefore explicitly retained both bosonic and fermionic 
degrees of freedom in our FRG calculation. The FRG 
approach developed in this work should also be useful to 
discuss other model systems where gapless fermionic and 
bosonic excitations are strongly coupled. 
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APPENDIX A: SKELETON EQUATIONS 

In sections |IV| and |VJ we have combined FRG flow 
equations with skeleton equations for the bosonic self- 
energy and the bosonic three-legged vertex to obtain a 
closed system of equations. In this appendix, we briefly 
describe the derivation of these skeleton equations. 

Skeleton equations relating vertex functions of differ- 
ent order follow from the general Dyson-Schwinger equa- 
tion for the generating functional for the connected Green 
function, which is a simple consequence of the invariance 
of the integration measure of the functional integral un- 
der infinitesimal shifts. The derivation of the skeleton 
equation for the bosonic self-energy of models of the type 
considered in this work has been discussed in detail in 
Refs. [36l and l43| so let us here only quote the result, 



■J K „, „ 



X Y%Y'^°''^"'^(K; K^K- -K). 



(Al) 



This exact identity is shown diagrammatically in Fig. [91 
To derive the skeleton equation for the three-boson ver- 
tex, let us start from the Dyson-Schwinger equation given 
in Eq. (11.27a) of Ref. [5S1 After taking two successive 
derivatives with respect to 



'K, 



and 



'Ks^ 



we obtain 



5^T 



Hk, h 



K2 



"K, 



ErS 



5^gc 



(A2) 
Here, T[tp,ip,(j)] is the generating functional of the one- 
line irreducible vertices, and Gc[j,j,J] is the generating 
functional of the connected Green functions, and is a 
functional of the sources j, j and J conjugate to the fields 
'0,'0, and (j). We now use Eq. (6.82) of Ref. l36land set 
all external fields equal to zero. The desired skeleton 
equation can then be written as 



T^^^{Ki,K2,-Ki -K2) = Ny2T^ f G''{K)G"{K + Ki)T^"^"'^'^{K; K + Ki;K2,-Ki - K2) 

Uk 

G°'{K)G"{K + Ki)G°'{K + Ki + K2) 

K 

^pV'°»A°0(X + Ki + K2; K + Ki;K2)T'l''''^'''''(K; K + Ki + K2; -K^ - K2) 
G"iK)G°'{K + Ki)G"{K - i^2)r'^°'^°'^(^ -K2;K + Ki; ~Ki ~ K2)V'^"'^'^'^(,K; K - K2; K2) 

K 



(A3) 
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y^ 




+ 




FIG. 10. (Color online) Graph of the skeleton equation (|A3[) for the three-legged bosonic vertex 



A graphical representation of this equation is given in Fig. [TO] Let us check a known limit of this equation. After 
replacing the exact propagators G"'{K) by Gq{K), the exact vertices V^ ^ '^{Ki; K2 ; K) by the bare vertices Fq , and 
neglecting the contribution involving F'^ '^ '^'^(Ki; K2:, Ki, K2), we obtain from Eq. (A3) 



Jk 



+ G^{K)G^iK + Ki)G^{K-K2) 
Jk 



(A4) 



which is invariant under arbitrary permutations of Ki, 
K2, and K3 = —Ki — K2. In fact, in this approximation, 
the bosonic three-legged vertex can be identified with the 
symmetrized closed fermion loop with three external legs 
and bare propagators given in Eq. (Bl|. As concerns 



Eq. ( A3 1 , the symmetry under permutations of its argu- 
ments is less obvious. However, as the left-hand side of 



Eq. (A2| is symmetric under permutations, and because 



all manipulations are exact, Eq. (A3 1 indeed fulfills this 
symmetry. 

To obtain the approximation for the bosonic three-loop 
used in Sec. Ivj we adopt the same approximation strat- 
egy as in the skeleton approximation used for the bosonic 
self-energy: replacing the vertices r'^'''^'''^{Ki]K2;K) 
and T'l''''f''''f"f'{Ki;K2;Ki,K2) by the bare vertices F^ 
and 0, but r etaini ng dr essed propagators G"{K), we ar- 
rive at Eqs. ( |5.13[ ) and ( |5.14[ ). 



Lf(iri,i?2,i?3), defined by 

L^{K,,K2,K3)^ f G^{K-K^)G^{K-K2)G^{K-K^) 
Jk 

(B2) 

as follows, 

L^iK,,K2,K3) 
= ^[L^iK,,K, + K2,0) + L^{K2, Ki + K2,0) 
+L^{K2, K2 + K3, 0) -I- L^iKs, K2 + X3, 0) 
-f Z^ (/?3 , i^3 + i^i , 0) -t- Z^ (Xi , i^3 + i^i , 0)] . (B3) 



Actually, taking into account energy-momentum conser- 
vation, we may set if 3 = —Ki — K2, so that we need 



APPENDIX B: THE SYMMETRIZED 
THREE-LOOP 



In the momentum-transfer cutoff scheme, all irre- 
ducible vertices involving only bosonic external legs are 
finite at the initial scale and can be identified by the sym- 
metrized closed fermion loops with bare fermionic prop- 
agators. In particular, the initial condition for the vertex 
with three external boson legs is 

r^**(i?i,i?2,i^3) = 2!iV^(F^)3i^(-/?l,-i?2,-i?3), 

a 

(Bl) 
where the symmetrized three-loop L"{Ki,K2,K3) can 
be expressed in terms of the non-symmetrized three-loop 



L^iK,,K2,-K,-K2) 

^[L^{Ki,K, + K2,0) + m^K,,-Ki - K2,0) 
+mK2,K, + K2,0) + m-K2,~K, - K2,0) 

+mKi,-K2,0) + L^i-Ki,K2,0)\. (B4) 



For our model, the three-loops can be calculated ana- 
lytically using the method outlined in the appendix of 
Ref. I3H1 Consider first the non-symmetrized three-loop 
defined in Eq. (B2). To perform the loop integration. 



we decompose the integrand in partial fractions, then 
carry out the ku integration by means of the method of 
residues, and finally perform the w-integration. Using 
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the notation Ki — {itOi, fc^), the result can be written as with 



'-^ij — l-Ui — Wj , Qij — kj_i — k±j 



ij(A-.,*.A,)=^/|^i:...n™i 



h — 



C_Li + Kj_j 



a{kni — kuj) 



^q^ 



(B7a) 
(B7b) 



1 r dkj^ 
2n / ^tT 



lUJi 



i=l 3 = 1 ^3 

iW2 



(fc±) 



The remaining fcj^-integration in Eq. ( B5 ) can now be 



' n^Ak±)^Uk±) ' ^ ^ 



done using the residue theorem and we finally obtain 



1 

'8^ 



where we have defined 



c^i e(Im(fci2)) - e(Im(fci3)) 
<7l29l3 ki2 - fci3 

Q2 e(Im(fc23)) ~ e(Im(fc2i)) 

923921 ^23 - ^21 

c^3 e(Im(fc3i))-e(Im(fc32)) " 

931932 ^31 - ^32 



.(B8) 



n?^ik^)^^u,~^u^,+^^_-^^-^^_-^^ 

= iujij - Q:(fc||i - fc||j) + {k]_^ - kj_j) - 2qijk± 

= ^Qij [n;± kij \ , 



Substituting this expression into Eq. ( B4 ) and defining 

k±i 



k±i 



(B6) 



k±i + k±2 ' 
we obtain for the symmetrized three-loop, 



(B9) 



L^iK,,K2,-K,-K2 



3\2tt 



S1S2 



feiii 



Ll K±2 



ia[xi — X2) 



k±ikj_2 



{si\xi\ +S2|x2|)e(-a;ia;2) 



-|-(si|xi| - \sixi + S2a;2|)0(-a;i(siXi + S2X2)) + {s2\x2\ ~ \siXi + S2a;2|)6(-a::2(sia;i + 522^2)) 



(BIO) 



Alternatively, this expression can be written as 
L^{Ki,K2,-Ki-K2 ^ 



3\2tt 



S1S2 



t-t-M-^--^ 



k±ik±2 



(siXi - S2X2)[Q{xi) - Q{X2)] 



-(2sia;i + 522^2) [0(a;i) - 6(siXi + S2X2)] + (2s2X2 + sia;i)[8(a;2) - 6(sia;i + S2X2)] 
1 siXi6(a;i) + S2a;26(a;2) - (sixi + S2X2)&isiXi + S2X2) 



47r 



S1S2 



k±ik±2 



TT s ia.\x\ — X2) 

^ /^ ^ ^^ -10 (fe) + -^e (i^) - (^1 + ^2)6 (^ 

47r VA:_Li /c_l2 






fell fc I 2 ^ fc I 



U}2 ' 
fe I 1 fe^2 ' 



(Bll) 



(fc±i + kx_2y 



For a different effective model for the nematic quan- 
tum critical point involving a quadratic energy dispersion 
(and hence a compact Fermi surface) the scaling prop- 
erties of fermion loops have recently been analyzed by 
Thier and MetzneiH®!, who found that also in this case the 
fermion loops exhibit a singular dependence on momenta 



and frequencies. However, to obtain consistent scaling 
properties of the effective interactions between bosonic 
fluctuations described by the fermion loops, one should 
use one-loop renormalized fermion propagators G^iK^ — 
iuj-z^" ^^ ^^^ loop integrations, which can be formally 
justified from the skeleton equation for the irreducible 
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vertex with three external bosonic legs, as discussed in 
Appendix A. The corresponding expression for the renor- 
malized three-loop can be obtained from Eq. (Bll) by 



replacing all external frequencies by w — >■ w/Z and mul- 
tiplying the loop by an overall factor of Z. The result is 
given in Eq. (|5.14|) . 
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